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ABSTRACT 


Solar Active Region NOAA 11158 has hosted a number of strong flares, in¬ 
cluding one X2.2 event. The current density and current helicity complexity 
properties are studied through cancellation analysis of their sign-singular mea¬ 
sure, which features power-law scaling. Spectral analysis is also performed, rev- 
elaing the presence of two separate scaling ranges with diffrent spectral index. 
The time evolution of parameters is discussed. Sudden changes of the cancella¬ 
tion exponents at the time of large flares, and the presence of correlation with 
EUV and X-ray flux, suggest that eruption of large flares can be linked to the 
small scale properties of the current structures. 


Subject headings: solar flares, active regions, solar physics 
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Introduction 


Solar magnetic 


such as solar flares (Benz 


activity is often accompanied by spectacular 


20081) and coronal mass ejections (Chen 2011). These eruptive, 


abrupt phenomena, 


highly energetic features can produce variations in the Sun-Earth connection, resulting for 


exampl e in geomagnet i c storms and ot 


ler disturbances that can affect human activities and 


health (Sc hwe nn 


2006 


Pulkkinen 


2007 ). The increasing amount and the improved quality 


of solar observations, both from Earth and from space, has provided enormous advances in 
the understanding of the physical processes occurring in the solar regions associated with 
flares, i.e. the solar active regions (ARs). These are regions where emerging photospheric 
magnetic field concentrates in bi- or multi-polar structures, which may include the presence 
of pores and sunspots. Driven by convective motions of the external layers of the Sun, 
active regions entangled and twisted magnetic fields can store a considerable amount of 
non-potential magnetic energy. The rapid release of such energy, probably due to dissipation 


in magnetic reconne c tion, is thought to be the basic mechanism of solar flares (ICarmichael 


1 96 4 


Sturrock 


1968; 


Hirayama 


197 


Kopp and Pneuman 


19761 ) 


One important topic in solar physics is the identification of magnetic signature of the 
occurrence of flares within ARs. In recent years, conv incing evidences of the turbulent nat ure 


of the ARs magnetic fields dynamics have emerged (Abramenko and Yurchvshvn 


2010b). 


In such context, the complexity of the photospheric magnetic field structure arises from the 


strongly nonlinear, coupled interactions between the plasma flow and magneticjluc 


on different sc ales, and results in the superposition of correlated structures (Frisch 


Biskamp 


nations 


1995 


9941) . Typical power-law spectra of photospheric magnetic energy have been 


reported ( Zhang et al 


phenomenology ([K olmogoro v 


2014). with 


1941 


spectr al indices compatible with a Kolmogorov-type 


Frisch 


1995 ). Turbulent ph otosp heric magnetic fielc 


have also been described as intermittent and multifractal (Abramenk o and Yur ch vshvn 
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2010b, 


20081 ). which are typical features of turbulent plasmas. This implies the presence of 


a hierarchy of correlated fluctuations, which concentrates energy on localized, small scale 
structures, where enhanced dissipation occurs. 


In order to capture the dynamical properties of the intermittent structures, cancellation 


analvsis has been recently used with the aim of correlating the complexity 

r of solar magnetic 

fields in active regions to 

lie occurrence of flares ( 

Yurchvshvn et al. 

2000 

Abramenko et al. 

1998; 

Sorriso-Valvo et al. 

2004; 

Yurchvshvn et al. 

2014 

). Such analysis has highlighted 


the importance of sign singularities in the energy storage process that could lead to flare 
eruption. In this paper we use for the first time cancellation analysis to describe the fine 
time and space resolution dynamics of AR NOAA 11158. Corroborated by the study of 
other observables, this analysis confirms the presence of nontrivial correlations between the 
topological changes of magnetic structures and flaring activity. 


2. Signed measure and cancellation analysis 


Solar active regions are often ch aracter i zed by scale dependent formation of energetic 


and localized magnetic structures ({Abramenko et al 


19981 ). Because of their coherence, 


structures can be seen as smooth regions of the magnetic field, embedded in a highly 
fluctuating background. For zero-mean fields, they can be associated to scale dependent, 
signed fluctuations of the fields. By introducing a signed measure (as opposed to the usual 
positive defined probability measure), it is possible to ch aracter ize the scaling properties 
of sign oscillations (or sign persistence) of the fields (lOtt et ah 1992). Therefore, the 
presence and the topological characteristics of structures defined in sign can be studied. 
Signed measure has been successf ully used to describe the cancellation properties of 


magnetic dynamo (lOtt et al 


19921) . as well as the characteristics of current structure s 


in turbulent magnetohydrodynamic (MHD) (Sorriso-Valvo et al. 


2002 


Graham et al. 
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2005 

), Hall-MHD ( 

Martin et al. 

2013 

). and kinetic ( 

De Vita et al. 

2014; 

Wan et al. 

2012) 


numerical simulations. Applications to measurements of magnetic vectors in solar active 
regions have confirmed cancellation analysis as an interesting tool to detect changes in the 


scaling pro pert ies of the fields flue 


gradients flYurchyshyn et ah 


2000 


uations, and of t 


Abramenko et ah 


re fractal dimension of the associated 


1998 


Sorriso-Valvo et ah 


2004). 


The signed measure of a zero-mean scalar held /(r) can be defined on a d -dimensional 


set Q(L) of size L as follows (Ott et al. 


1992). Let { Qi(l ) C Q{L)} be a partition of Q(L ) 


in disjoint subsets of size l. Then, for each scale l and for each disjoint set of boxes Qi(l), 
the signed measure is 


/h(0 = 


4(0 dr ^ r ) 


( 1 ) 


Iq(l) dr l/WI ' 

When the size of the subset Qi(l) is large, cancellations between small structures of opposite 
sign occur within each box, resulting in small contribution to the signed measure. However, 
as the boxes become smaller and reach the typical size of the structures, each box is more 
likely to contain one single, sign defined structure, reducing the level of cancellations. The 
way this happens can be statistically characterized through the cancellation function 


X (i) = £ mq i 

Qi(l) 


( 2 ) 


where the sum is extended to all disjoint subsets Qi(l). Contrary to positive defined 
probability measures, the signed measure holds information on the sign of the held 
fluctuations. In particular, if the measure changes sign on arbitrarily fine scale (i.e. if for 
any subset Qa(1) for which 7^ 0 there exists a subset Qb( l') C Qa(1) such that ur(I') 
has opposite sign from ha(1)), then the measure is called sign-singular (IQtt et al.lIl992h . 


Upon performing a scale dependent partition of the whole domain, the sign-singularity of 
the measure can be quantitatively estimated through the cancellation exponent k, that is 
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the scaling exponent of the cancellation function, defined as 


x(0 = £ MOI ~ i- K , 

Qi(l ) 


( 3 ) 


where the sum is extended to all disjoint subsets Qi(l) that cover the domain Q(L). In 
a fluctuating field with positive and negative structures, cancellations occur in large size 
subsets, providing small contribution to the signed measure. Conversely, when the subsets 
become smaller and reach the typical size of the structures, the enhanced presence of 
sign-defined structures reduces the level of cancellations, and consequently increase the 
relative contribution to the signed measure. Thus, the cancellation exponent represents an 
effective measure of the efficiency of the field cancellations. Specific examples are represented 
by a smooth field (with no sign-singularity), for which the cancellation function has a 
constant value (so that k — 0), and by a homogeneous held with random discontinuities (i.e. 
Brownian noise), for which n = d/2. Cancellation exponents between those two limiting 
values indicate presence of smooth structures embedded in random fluctuations. Moreover, 
their values can be related to the geometrical properties of structures. Furthermore, it is 


possible to establish a phenomenological relationship between the cancel 
the fractal dimension D of the typical dissipative structures of the how 


ation exponent and 


Sorriso-Valvo et al. 


20 021 ). 


n = (d — D )/2 . 


( 4 ) 


On the other hand, cancellation exponents larger than the typical value for random 
helds k > d/2 indicate the presence of sign anti-persistent structures. Such behaviour is 
typically associated to the presence of pairs of adjacent, opposite sign structures, which 
enhance cancellations with respect to random fluctuations. It should be pointed out that 
the use of one single fractal dimension cannot fully capture all the fine details of the 
typical plasma tubulcnce pr ocesse s, which are more likely characterized by multifractal 


scaling ( Muller and Bi skam p 2000). Nonetheless, D still represents a useful indicator of the 
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topological characteristics of the “mean” intermittent structures of the flow. 


3. Solar Active Region NOAA 11158 
3.1. Data reduction: AR 11158 


We derive evolution of magnetic fields in NOAA 11158 using series of HMI vector 
magnetograms. A 6-day uninterrupted, 12-minute cadence data-set allowed us to study in 
detail both the long-term, gradual evolution, as well as the rapid changes during an X-class 
flare. In this section we describe the data-set and the re-projection we use. 


NOAA 11158 was the source of an X2.2 flare on 2011/02/15 starting at 01:44 UT, 
peakin g at 01:56 UT and ending at 02:06 UT. A front-side halo CME accompanied the 


flare (jSchriiver et al. 


201 1). Prior to the X2.2 flare, the largest flare in this region was an 


M6.6 on 2011/02/13 at 17:28 UT, a little more than 30 hours before the flare of study. 

Helioseismic and Magnet ic Imager (HMI) of NASA’s Solar Dynamics Observatory 
(SDO) satellite ( Pesnell et ah 2012) observed the NOAA 11158 in high detail, routinely 
generating hltergrams in six polarization states at six wavelengths on the Fe I 617.3 nm 
spectral line. From these hltergrams, images for the Stokes parameters, I, Q, U, 
and V were deri ved w hich, using the Very Fast Inversion of the Stokes Algorithm 


(VFISV) code (Borrero et al 
components. To resolve the 


2 0111 ). were inverted into the magnetic held vector 


method ( I Met calf et ah 


1994; 


80° azirnut 


Leka et ah 


ral held ambiguity we used the “minimum energy” 


2009 ). In addition, we corrected a few episodes of 


single-frame fluctuations in the direction of the transverse magnetic held vector, by nearly 


180° ( Welsch et al. 


2013 1). 


To study pre-hare photospheric magnetic evolution, and to baseline this evolution 


against post-hare evolution, we retained 153 hours of 12-minute-cadence 0.5”-pixel 


























resolution HMI vector-magnetogram data, from beginning of the active region emergence, 


around four days before the X2.2 flare, to two days after the flare: t start = February 10 2011 
14:11 UT (519,^50), t end = February 16 2011 23:35 UT (S21,W37). 


To account for the non-zero viewing angle of the observed vector magnetic fields and 


Doppler velocities, we repro jected t he 


to Cartesian coordinates flWelsch et al. 


Mi-data cube to the disk center and transformed it 


20131) . To do that, in the first step, we re-projected 


the observed magnetic vectors’ components onto radial/horizontal coordinate axes. We 
then determined the shift of the grid’s center between frames. After we determined a shift, 
we converted and then translated by the decided-npon shift the Cartesian grid’s points 
to LOS/POS coordinates. We then interpolated the radial and horizontal components 
of the magnetic field, B r , B h , onto the grid. Finally, we interpolated data in LOS/POS 
coords onto a fixed, Cartesian grid that accounts for shifts. This interpolation implied a 
re-projection of the magnetogram surface, from spherical to a regular, Cartesian grid. Since 
we used Fourier Local Correlation Tracking (FLCT) to find velocities, we ne eded co nformal 


trans formation into cartesian coordinates; we chose Mercator re-projection i lWelsch et al. 


20091 ). with regular longitudes/latitudes. After re-projection, to preserve physical quantities 


of magnetic fields and velocities, we corrected the fluxes for the distortion of pixel areas 


introduced by re-projection; the c 


Appendix A of 


Kazachenko et al. 


etails of the applied correction-factors are given in 


2014bh . For the minimum magnetic field to consider, we 


chose a threshold of |B| = 250 G auss, c onsistent with the upper limit of the uncertainty in 


the horizontal magnetic field (jLi u et al. 


20121 ). To avoid spurious signals in magnetic fields, 


we apply a mask, where we set any pixel’s magnetic field components to zero, if in any of 
the three consecutive frames it has |B| < 250 G auss. We also added a boundary area of 


55-pixels width/height of padded with zeroes (Kazachenko et al, 20T4aj). The final data 


cube after re-projection and boundary adding consists of 768 time steps (dt = 720s) and 
has a field of view of 665 x 645 pixels with a pixel size of 360.16 km, which is equivalent to 
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the original 0.5” pixels HMI-data resolution. More details on th e data c ub e p reparation 


and calibration, for a shorter time range, could be found in 


Welsch et al 


(2013). 


Fig. |T] shows the final LOS magnetic field in a subregion of the full-disk data array 
after re-projection in the beginning (Panel A), middle ( B,C ) and the post-flare (. D ) times 
of the magnetogram sequence. Note that the positive and negative magnetic fluxes, shown 
on the right plot, balance each other; the magnetic flux increases from essentially zero to 
roughly 1.4 x 10 22 Mx at the time of the flare (vertical dashed line). To summarize, as a 


B.: 201 1.02.10 14:1 1 B.: 2011.02.12 01:59 B.: 201 1.02.13 1 3:59 B,: 201 1.02.1 5 01:59 



z Pos/Neg Magnetic 
2000 fluxes, [iCPMxl 



- 2 ooo 11 12 13 14 15 16 
February 2011, [days] 


Fig. 1.— Panels A-D: HMI vertical magnetic field ( B z ) maps at 4 different times of AR 
11158 evolution. Panel E : evolution of the positive and negative vertical magnetic fluxes of 
the 6-day interval, with the diamonds indicating the times of the four images (A-D) on the 
left. An X-class flare occurred at the time corresponding to the vertical dashed line. The 
black box in Panel D indicates the Region where the X2.2 flare was located. 

result of the data reduction we obtained a data cube consisting of 768 frames, each of which 
contains data for three components of the magnetic held with a held of view of 665 x 645 
pixels, a pixel size of 360.16 km and a time step of dt = 720 sec. 

For sake of complcnteness, Figure [2] shows the three magnetic held components at 
two different times, on February 12th at 22:24UT (left panels) and on February 14th at 
03:24UT (right panels), before and after the hux emergence. The black solid line box 
indicates the reduced area used for the cancellation analysis. From the hgure, the evolution 
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of the magnetic structure of the AR is evident. The full time evolution of the magnetic 
field component B z is shown in the movie 13.11 where it is possible to observe the clear flux 
emergence, and the successive increase of the magnetic complexity. 

2011 . Feb 12 , 22:24 UT 2011 , Feb 14 , 03:24 UT 



Fig. 2.— The magnetic field components intensity measured for NOAA 11158 on February 
12th at 22:24UT (left panels) and on February 14th at 03:24UT (right panels). 


I\ Gauss 

























































11 


Fig. 3.— A clip of the time evolution of the dataset and relative quantities evaluated in this 
paper. 

3.2. Preliminary analysis and complementary measurements 


The measurements of the magnetic field vector allow to calculate the current 
density component perpendicular to the solar surface, J z , shown in panels (a) and (b) 
of Figure [4] at the two times presented in Figure [2l The vertical component of the 
current density J z (x,y) has been estimated through the photospheric vector magnetic 
field B (x,y), where (x,y) are the cartesian coordinates on the solar surface, as the 
line integral of the transversal component of the magnetic field over a closed contour 
G, AnJ-jc = (V x B). = s-'&Bx • dr. The integral along each pixel G of area 
s = (0.361 x 0.361 )Mm 2 was computed using Simpson’s rule. An alternative calculation, 
based on the finite differences evaluation of the magnetic field rotational, has been also 
performed, providing identical results. Previous studies of cancellations in solar AR have 


shown that it is conve nient to use the reduced current helicity h c = J Z B Z (lAbramenko et al. 


1998; 


Jing et al. 


2012). This field carries important information about the non-potential 


magnetic energy available in the AR. Furthermore, it is normally less noisy than the current 
density, so that cancellation effects are easier to identify. As a matter of fact, previous 
analyses were unable to identify sign-singularities in the current density, which were only 
observed in current helicity. It should be pointed out that, since it is not possible to measure 
the full current density vector with the available magnetic field measurements, in this work 
we limit our discussion to the current vertical component and the reduced current helicity 
h c (x,y) = B z (x,y)J z (x,y), and assume that they are approximating the corresponding full 
quantities. Panels (c) and (d) of figure Q] show, for the same two snapshots as in previous 
figures, the estimated current helicity, which has a smoother appearence than the current 
J z , as expected. Movie I3TT1 reproduces the full temporal evolution of the current helicity 
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field, highlighting the important establishment of the AR current structure at one given 
time (see below). 


a) J z 








Fig. 4.— The vertical component of the current density J z (panels a, b) and the reduced 
current helicity h c (panels c, d), computed for NOAA 11158 on February 12th at 22:24UT 
(left panels) and on February 14th at 03:24UT (right panels). 

The amount of power connected to small scale magnetic fields can be quantitatively 
visualized through the typical fluctuation level of the magnetic field vector, B rms = 
(AB 2 ) 1//2 + ( 8B 2 ) 1//2 + (hR 2 ) 1 / 2 , where SBi = B\ — ( B { ) are the fluctuations of the i-th 
magnetic field component, and the brakets indicate space average over the whole AR. 
Figure [5b shows the temporal evolution of B rms , indicating the increase of fluctuations 
starting at the beginning of the emerging phase of the AR (see the magnetic flux), and then 
an evident decrease in the final part of the obsevation. 


fiA/rri' 
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Another useful quantity is the mean squared current (Jf), the average being again 
estimated over either the whole AR. This quantity is related to the level of dissipation of 
magnetic energy, and is customarily used in numerical simulations to describe the evolution 
of turbulence. Panel (c) of Figure [5] shows the time profile of (J|) (black solid line). From 
the figure, a sharp increase of (Jf) is seen at t* =1:24 UT on February 13th, during the 
initial stage of the flux emergence. This interestingly suggests that, unlike the smooth, 
slower and progressive emergence of magnetic flux and magnetic fluctuations, the transition 
toward a high mean vertical current state occurs in a very short time, At* = 192 minutes. 

Finally, since we aim at connecting the magnetic properties of the AR with occurrence 
of flares, we have collected the measurements of the X-ray (X) and Extreme-Ultra-Violet 
(EUV) fluxes, as measured by GOES and AIA/SDO (13lA channel), respectively 
(Figure EK). 

GOES provides the X-ray flux measurements integrated over the whole solar disc (black 
solid line). This could include in the measurement flaring activity from other active regions. 
However, we have checked that during the time interval under study (the gray shaded area 
in all panels of Figure [5]), the X-ray flux was mainly originated in NOAA 11158. 

On the contrary, the EUV flux measurement retains the spatial information. Images 
taken at the 13lA band pass have been selected for their sensitivity to hot (~ 10MK) 
plasma to complement the GOES data. Level 1.5 data (cutout service) have been used at 
a 12-minute cadence. The maps have been integrated over the active region NOAA 11158, 
resulting in a time series of EUV emission from the area under study (red dot-dashed 
line in Figure [5]). Some data gaps are present in the EUV time series. For the statistical 
analysis, such gaps have been filled through cubic spline interpolation of the data. During 
flaring time, the brightest parts of the flaring region are generally saturated. Therefore, the 
absolute value of the EUV emission does not correctly represents the intensity of the flare. 
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Nevertheless, higher intensity might still represent larger flares, as generally large flares 
tend to have larger areas of saturated pixels. In any case, what is important for the study 
present in this paper is that the time series indicates flaring times properly. Therefore the 
saturated pixels are not of importance for the conclusions drawn in this paper. Moreover, 
the correspondence with flares eruptions estimated through the X-ray flux avoids any 
ambiguity. 

Figure |5^, shows that, after a quiet interval in the early stage of the AR prior to 
the flux emergence, both X and EUV flux background level sensibly increases, starting 
approximately at the transition time t*, in agrement with the increase of the dissipation 
proxy ( J 2 ). After this time, several flares were released by NOAA 11158, as evident from 
the peaks in the temporal profile of both X-ray and EUV flux. These include 18 C-class, 
three M-class and one X2.2-class flares. Vertical dashed lines identify the exact time of 
eruption of the M and X flares in all panels of Figure [5] 


4. Scaling laws in NOAA 11158: spectral and cancellation analysis 


4.1. Spectral analysis 


In order to characterize the possible presence of turbulent-like behaviours in the 
active region magnetic field, as suggested by the large mean current density, for each 
snapshot of the time series we have calculated the (omnidirectional) ma gne tic spectrum 


E B (k ) = Jj k | =A . |R(k)| 2 dk, where k indicates the wavevector fjZhang et all 120141) . Examples 


/|k| 

of magnetic spectra are shown in Figure El at the same two times as in Figured suggesting 
the presence of a power-law decay in the intermediate range of scales (the inertial range), 
approximately between I. 8 M 111 and 6 Mm. In this range, a power-law fit Ef,(k) oc k~ aiar s e 
provides the scale spectral index ai arge . At the early stage of the AR emersion, the spectral 
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power is smaller (Figure EJ), the power-law is less defined and the spectral index is very 
variable. At later times t > t*, when the active region is emerging, the scaling exponent is 
more stable, a ~ 2. This behaviour is clearly visible in the supplemental material 13.11 where 
the movie reproduces the time evolution of the spectra. The observation of a power-law 
spectrum, with spectral index compatible with a Kolmogorov-like phenomenology, confirms 


once again that the AR magnetic fields can be studied in the 


already suggested in the past (Abramenko and Yurchvshv: 


ramew ork of turbulent flows, as 


2010b). At larger wavevectors, 


the spectrum is compatible with the presence of a secondary, different power-law, with 
a sma u ~ 3.3, although the range of scales is rather limited. This behaviour indicates the 
presence of a characteristic scale, around 1.8Mm, where the physical processes change. Note 
that recent esti maties of the Batchelor integral scales in the quiet solar photosphere provided 


similar values ({Abramenko 


2013 1. The presence of a spectral break is commonly observed 


in plasmas, for example at the transition from the magnetohydro dynamic range and the 
kinetic scales, where particle effects are not negligible and produces variations in the energy 


cascading mechanisms, resulting in steeper spectra. 


Examples of spectra 


steeper sp ectra are o 


2007; 


20071: 


Sahraoui et al 


Wanetal 


rserve d in solar wind plasmas (ILeamon et al 


1998 


break followed by 


Alexandrova et al 


2009 ). as well as in a variety of numerical simulations (Mininni et al 


2012 ). At very large wave vectors the spectrum flattens, suggesting that 


scales < IMm may be affected by instrumental noise. These s pectr al propertie s are ro ughly 


in ag r ement with previous analysis of the sam e acti ve region (IGeorgoul is 


2014 


Guerra et al. 


2013; 


2014) and of other ARs (Abramenko and Yurchvshvn 


Zhang et al. 


2010 ai ) where, 


however, only inertial range spectal indices were evaluated. 
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4.2. Cancellation analysis 


Upon confirmation that the active region magnetic field can be described in the 
framework of turbulence, cancellation analysis has been performed on the vertical current 
J z and on the reduced current hclicity h c , for each snapshot of the time series. Note that 


the magnetic field is relatively smooth, so that i ts cance llation ana. 


information about the dynamics of the AR ([Abramenko et ah 


ysis does not provide 


1 998 ). Figure [7] shows one 


example of the signed measure maps, calculated from the data through equation (JTJ) at 
four different partition scales. While at large partition scales the positive and negative 
fluctuations cancel each other, resulting in small values of the signed measure, at smaller 
partition scales the sign-defined structures emerge as brighter regions of the maps. Note 
also the sign-defined large structures in the umbral part of the sunspots, where the magnetic 
field direction is well-defined. The way structures influence the overall signed measure 
at different scales is resumed by the cancellation function, depicted (for the same two 
snapshots as in previus figures) in Figure [HI for both the current (top panels) and the 
current helicity (bottom panels). A clear power-law range is visible in the intermediate 
range of scales, conservatively between 1.4Mm and 6Mm, roughly corresponding to the 
“inertial range” of the spectrum shown in Figure [6] A power-law fit has been performed in 
such range, as indicated in Figure [HI and the corresponding cancellation exponents have 
been evaluated. For the examples given here, exponents are steeper before the transition 
time, Kj z = 0.72 ± 0.02 and Kh c = 0.55 ± 0.03, while after the transition they become 
shallower, kj z = 0.41 ± 0.01 and Kh c = 0.13 ± 0.01. Based on the phenomenological 
argument given in Section [2j these values correspond to the presence, in the emerging 
stage, of broken current filaments with fractal dimension Dj z = 0.56 ± 0.04, and current 
hclicity filaments, Dh c = 0.9 ± 0.06. At later stage, current filaments dominate, with fractal 
dimension Dj z = 1.18 ± 0.02, while well resolved, almost smooth current hclicity structures 


are observed, Dh c = 1.74 ± 0.02. The faint saturation to y = 1 at small scales shows 
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that small scale features are rather smooth, probably becuase of the instrumental noise, 
and in agrement with the spectral observations. The extremely high quality of HMI-SDO 
data thus allows a very accurate estimation of the cancellation effect. The complete 
temporal evolution of the cancellation function is represented in the movie 13.11 where it is 
possible to confirm the good quality of the fit during the whole AR observation, and the 
clear changes occurring during the AR evolution. As a stricking difference with previous 
results, sign singularities are measured not only in the smoother current helicity, but also 
directly on the current density. This is the first observation of this kind, witnessing the 
excellent quality of the data, and providing further information on the fields dynamics. 


Note also that, w’ 


rile earlyer results on AR cancellation ana 


space resolution (Abr amenk o et al. 


1998 


Yurchyshvn et al. 


ysis were affected by limited 


2000 


Sorriso-Valvo et al. 


20041 ). a more recent study o: 


resolution flYurchyshyn et al. 
resolution. 


Hinode magnetic fields has shown a similar excellent spatial 


2014), but was still lacking full time coverage at high temporal 


5. Results: time evolution of the magnetic field complexity and relationship 

with large flares 


Previous results on cancellation analysis have suggested that the topological properties 


of active regions magnetic fielc 

s experience changes in corres 

>pondence to eruption of 


ma i or flares ( 

Abramenko et al. 

1998; 

Yurchvshvn et al. 

2000; 

Sorriso-Valvo et al. 

2004 


Yurchvshvn et al. 

2014 

). The first observations indicated that such abrupt changes may 


anticipate the flares by a fraction of a hour. However, the limited time resolution of 
the data used so far has been a serious obstacle to the detailed description of the time 
evolution of the magnetic field complexity. HMI-SDO data used here provide for the 
first time high time resolution, and full time coverage during almost the whole lifetime of 
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the active region. This dataset represents a unique opportunity to study the dynamical 
properties of the magnetic complexity, and their relationship with occurrence of flares. 
Previous studies of temporal evolution of AR NOAA 11158 focused on the properties of 


magnetic flux, potential, non potential and free energy, num 


Der o 


Tziotziou et al. 

2013; 

Liu and Schuck 

2012) 


Aschwanden et ah 


,2014 


l oops, m agnet ic hclicity, 


Vemareddv et al. 


2012 


5.1. General temporal dynamical properties: turbulent dissipation and flares 

Looking at the gross behaviour of the parameters can give information about the 
general magnetic and energetic properties of the active region. To this aim, in Figure [5] 
we have collected the time-dependent parameters of the active region, as measured by the 
instruments or calculated for this work. We have already discussed in previous sections the 
general properties of the time evolution of: (a) the EUV and X-ray flux, useful to identify 
the time of eruption of flares calssihcation (the time of eruption of three M and one X flares 
is indicated by vertical dot-dashed black lines); (b) the positive ( F pos ) and negative (F neg ) 
magnetic fluxes of the whole active region, indicating the slow and steady emergence of 
magnetic flux, and the similar behaviour of the magnetic fluctuations level ( B rms ); (c) the 
space-averaged squared vertical current (Jf), which shows the abrupt onset of turbulence 
at time t*, and a successive steady state (indicated as a gray area in the plot) preceding a 
final smooth decrease. The other panels of Figure [5] show: (d) the spectral indices a sma ii 
(thin blue line) and cq arffe (thick yellow line), obtained from the power-law fit of the spectra 
as described in Section [4] for both ranges of scale; (e) the cancellation exponents estimated 
from the power-law fit of the cancellation functions, for both the vertical current (ftj z , green 
thin line) and the reduced current helicity (nh c , thick red line). The complete temporal 
evolution of the cancellation exponents can be observed in more detail in the movie 13.11 
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The time evolution of the spectral indices, shown in panel (d), is only weakly affected 
by the AR dynamics. Indeed, the spectral exponent reaches a steady state at early times 
(about February 12th for the small range spectral index, and mid February 11th for the 
large scale exponent), one-two days in advance with respect to the flux emergence and the 
sharp jump observed in (Jf). While there is no evidence of later evolution for ai arg e (for 


example at the f ull emersion of t 


flaring activity) flGeorgoulis 


re_AR, or in cor resp ondence with the enhancement of the 


2013 


Guerra et al. 


2014), the small scale spectral index seems 


to start increasing at t*, and then reaches a broad steady state, starting from the beginning 
of day 14th up to about 14 hours in day 15th (gray area in Figure [^D- The larger values 
a sma ii ~ 3.5 indicate a steeper scaling exponent, which could be attributed to a time interval 
of more efficient transport and dissipation of turbulent energy. This spectral modification, 
to our knowledge observed here for the first time, is in agreement with the increase of X-ray 
and EUV background flux, as well as with the enhanced flaring activity. In either cases, 
while the large scale properties of the AR magnetic fluctuations do not seem to be affected 
by the flaring activity, the small scales are sensitive (although just weakly) to the erupting 
phase of the AR. This observation confirms that flares are most likely connected to small 
scale dynamics. 

Finally, the time evolution of the cancellation exponents is shown in panel (e). In 
the early stage of the AR emergence, the current and current helicity structures are not 
well defined yet, the AR being dominated by noise. In these conditions, the cancellation 
functions show variable sign-singularity, so that the cancellation exponents are highly 
fluctuating. A similar behaviour was recently observed in direct numerical simulations for 


the s t udy of the trans ition to turbulence in kinetic dynamics of plasmas (IDe Vita et al. 


2014 


Wan et al. 


2012J), and is due to the presence of not yet fully developed structures. 


At t*, an abrupt decrease of the exponents synchronizes well with the sharp jump of the 
averaged squared current density. This indicates, once again, that the complexity of the 
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magnetic field changes in a very short time, according to the fast increase of currents in the 
AR. This change is followed by a slower, steady decrase, lasting about one day. Then, the 
two cancellation exponents reach a steady state, with values Kj z ~ 0.42 and hih c — 0.13, 
typical of smoother structures expected in a highly dissipating plasma (gray area in the 
plot). After about 36 hours, the cancellation exponents start growing again, indicating an 
increase of complexity of the AR magnetic fields. 

The whole gross temporal evolution of the AR can therefore be resumed as follows: (1) 
the eary stage of the AR, with no flaring activity, showing randomly emeging, disrupted 
current filaments associated with weak turbulent energy; (2) the setup of the flaring 
activity and magnetic flux emergence, marked by the sharp onset of turbulence occurring 
at t*, with stabilization of the structures, and the following day of steadier change in the 
paramenters; (3) a period of strong flaring activity, associated with enhanced level of 
magnetic fluctuations and turbulent dissipation, and charaterized by a steady state of the 
structures geometry (gray area in Figure |5]); (4) finally, the weakening of the flaring activity, 
associated with the steady decrease of turbulence level, weakening of the dissipation 
(shallower small scale spectrum), and increase of the magnetic complexity. This final step 
of the AR temporal evolution seems to indicate the transition to a different state, where 
finer magnetic structures can build up without necessarily result in large flares, suggesting 
an improved capacity of energy storage at smaller scales. The study of this part of the 
evolution is however not going to be discussed any further in this paper. This scenario is 
fully consistent with the association between eruption of flares, and the general properties 
of the AR dissipation and small scale magnetic complexity. 
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5.2. Short time-scale features at flaring times 


Besides the gross dynamical behaviour of the AR, the high temporal resolution of HMI- 
SDO allows for the first time the detailed analysis of short time-scale features, which are 
extremely suitable for the study of rapid phenomena such as solar flares. The instrumental 
cadence of 12 minutes is in fact comparable with the typical time scales associated to the 
flaring process. In order to study the possible relationship bet wee n f lares and m agnetic 


turbulence properties, already pointed out in the recent years flSorriso-Valvo et al. 


2004 ). 


the study can be focused on the most intense flares occurred in NOAA 11158. Previous 
studies have indeed shown that mostly major flares can be associated with variations of the 
magnetic complexity, as can be expected. In Figure 0 four vertical lines indicate the time 
of eruption of the largest flares of NOAA 11158. As already mentioned, these are three 
M-class (one of which considerably larger than the others) and one X-class flare. 

It is evident from the panels (b) and (d) of Figure 0 that, at the time of the flares, no 
specific features are observed in the magnetic flux and in the magnetic spectral exponents. 
However, finer observation of such and other similar parameters evalu ated in th e flare 


triggering regions have revealed small but significative changes (Petrie 


2013; 


Song et al. 


2013). On the contrary, the magnetic fluctuations B rms show a clear sudden increase at the 


times of the first M and X2.2 flares. This is highlighted in Figure [9j where a magnification 
around the M and X flaring periods is shown (panel a). In the same figure, the time 


derivative of the two fields are also shown, with the evident peaks at the first M and at the 

X flares. No features are visible for the second and fourth M flares ( 

Petrie 

2013; 

Sons: et al. 


2013). At the same times, the mean squared current (Jf) also exhibits interesting features. 


The M-class flare is associated to a jump in the mean dissipation, while at the X-class 
flare an evident broad peak is present. In the latter case, the turbulent dissipation level 
increases by 27% (two orders of magnitude larger than the 0.2% relative standard deviation 
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evaluated in the steady period preceding the flare) in about 216 minutes, peaks about 


48 minutes before t 


re flare, and then decreases back to its steady value in about the 


same amount of time ([Song et al. 


20131 ). Similar features are also observed at the X flare 


for both cancellation exponents, where a 6% and 28% growth is present for the current 
density and current helicity indicators, respectively. These increases are larger, by one 
order of magnitude, than the 0.5% and 0.8% standard deviation levels in the steady period 
preceding the flare. The duration of the increase and decrease phases is approximately the 
same as for (J|), i.e. of the order of 200 minutes. The same kind of behaviour seems to 
hold for the weaker features observed at the time of the first M-class flare for the current 
cancellation exponent, although these changes are not as evident as for the X2.2 flare, and 
are absent in the current helicity exponent. 


The increase of the mean squared current and of the magnetic complexity (as estimated 
through the cancellation exponents), which represent the main result of this work, can be 
interpreted as follows. During that phase, the AR is in a highly turbulent state (as shown 
by the spectra), with a steady, high level of dissipation (shown by the slightly steeper small 
scale spectra with respect to the non-flaring stage, and by the higher background of X-ray 
and EUV flux), and associated with the presence of relatively smooth current and current 
helicity structures. Shortly before the flare, magnetic gradients and complexity increase, 
suggesting the injection of an excess of magnetic energy, which is not fully dissipated, but 
rather stored through the build-up of field complexity, with resulting enhancement of the 
current filamentation in the AR. This finally results in the conditions for flaring. After 
the flare explosion, the conditions come back to the steady state, while the typical level of 
magnetic fluctuations starts to decrease. 


It should be mentioned that the cancellation exponent features shown here are opposite 
with respect to previous cancellation analysis, where complexity was observed to decrease 
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before a large flare (Abramenko et al. 


2004 


Yurchvshvn et al 


1998 


Yurchvshvn et al 


2000 ; 


Sorriso-Valvo et ah 


2014). This difference could be due both to the lack of good time 


resolution, which could prevent the correct temporal description of the parameters, and to 
the poor spatial resolution of previous analysis. 


6. Evaluating correlations 

The observation of the temporal behaviour of the active region magnetic structure has 
revealed the correspondence between changes in magnetic field properties, and the eruption 
of large flares. This is visible both in the gross evolution of the parameters and in the short 
time-scale features, as shown in Figure El In order to give a quantitative measure of the 
relationship between the observed features and the occurrence of flares, a statistical study 
of the correlations between the different parameters presented in Figure [5] can be performed. 
This is possible thanks to the high time cadence of the data, which allows for the first time 
a significant statistical study. In order to limit the effect of transients, noise and smaller 
flares, we performed the correlation analysis only in the stationary stage of the activity, 
indicated by the gray area in Figure El The cross-correlation coefficients between several 
pairs of parameters has been computed, and results are collected in Table El We have 
used both the Pearson and the Spearman (ranks) cross-correlation coefficients, evaluated 
at the time lag where they are larger (in Table El time lag of 0 hours indicate absence 
of significative correlation). The Spearman coefficient is less sensitive to nonstationarity 
of the samples, and to nonlinearity of correlations. Since there is no particular reason to 
expect linear correlations, we believe that Spearman coefficient could be a more efficent 
parameter accounting for correlations in this system. Examples of correlation functions 
(for both ordinary and Spearman ranks correlations) are shown in left panels (a, c, e) 
of Figure [10l Peaks are evident for the current helicity cancellation exponents and for 
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the small scale spectral index, while correlations are poor for the mean squared vertical 
current. The dashed vertical line indicates the time lag at which the Spearman correlation 
is maximum. As can be seen in Table [HI while some of the pairs are not correlated, for 
others it is possible to highlight a relevant correlation. For example, EUV flux exhibit a 
strong degree of correlation ps = 0.59 with the current helicity complexity Kh c (at r = —1.4 
hours lag), or p s = 0.55 with the small scale spectral index a sma u (at r = 1.8 hours lag). 
Note that maximum correlation time lags are negative, indicating that the changes in the 
parameters anticipate the X-ray and EUV emission, except for the spectral indices, which 
on the contrary react after the flares. The sign of the correlation suggests the causality 
direction of our observations. Increase of magnetic complexity and gradients anticipates 
the flares, while the dissipation becomes more efficient after the flares. The most relevant 
correlation is found between EUV flux and current elicity cancellation exponent. This 
confirms that sign-singularity analysis is a suitable, sensitive tool, able to capture the fine 
variations in the AR magnetic complexity preceding the eruption of large flares. 

Finally, more details about the nature of the observed correlations can be evidenced by 
showing the corresponding scatter plots for each pair of parameters, spaced with the time 
lag of the Spearman correlation peaks. These are shown in the right panels (b, d, f) of 
Figure (TUJ It appears evident that most of the correlation comes from the large flares (the 
top part of the scatter plots), while smaller X or EUV records (the background emission) 
are more randomly distributed with the other parameters. Interestingly, all major flares 
occur above given values of the correlated variables, i.e. for the cancellation exponents 
and for the small scale spectral index. This suggests the presence of a threshold of the 
magnetic dissipation and topological complexity, below which no large flares are observed. 
This quantitative observation further confirms that flares are strictly related to small scale 
turbulent and dissipative processes in the photospheric magnetic fields, and that higher 
complexity of the currents enhances the probability of observing large flares. 
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Table 1: Maximum correlation coefficients pp and ps, and the time lags Tp and Ts (in hours), 
estimated between the X-ray and EUV fluxes and the indicated parameters. 
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7. Conclusions 

Motivated by the extremely good quality of HMI-SDO photospheric magnetic field 
vector measurements, we have studied some magnetic properties of the solar active region 
NOAA 11158, and their relationship with eruption of large flares (M and X class). We have 
observed that the dynamics of magnetic fluctuations and mean vertical current describes 
quite well the transition of the AR into the flare activity stage. The spectral properties of 
the magnetic field fluctuations were also studied, suggesting the presence of a double range 
of scales. In the larger scales range, approximately corresponding to the inertial range 
of turbulence, the spectral properties are steady during most of the AR lifetime, and no 
particular features are observed at the time of major flares. On the contrary, the small 
scale range is characterized by a variable spectral index, which shows nontrivial correlations 
with the X-ray and EUV emission. In particular, the spectral slope increases (indicating 
more efficient dissipation of turbulent energy) during the central part of the observation, 
when flaring activity is enhanced. The magnetic field complexity was studied by means 
of the cancellation analysis of its sign-singularities, as evidenced through the study of the 
vertical current and of the reduced current helicity. Cancellation analysis applied to the 
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AR provided the qualitative estimation of the fractal dimension of the current structures, 
which is reasonably steady during the flaring time interval. However, interesting peaks are 
observed in correspondence of the largest flare of the AR (X2.2), and to some extent also 
for the M class flares. The results shown here suggests that it is possible to quantitatively 
measure the magnetic complexity evolution of the AR during flares, which support the 
scenario of increasing entanglement of current and magnetic held one-two hours ahead of 
big flares. 

Any comments on the consequences on the main flaring models? Should we 
mention that there is no need to focus at the location of the flare, as necessary 
for other observables (see literature)?. 


The research leading to these results has received funding from the European 
Community Seventh Framework Programme (FP7/2007-2013) under grant agreement n. 
269297/ “TURBOPLASMAS”. 
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Fig. 5.— Time evolution of: (a) flux of X-ray (black full line) and EUV (red dot-dashed 
line); (b) space-integrated positive (thin, full blue line) and negative (thick, red dashed line) 
magnetic flux, together with the level of magnetic fluctuations (thick solid black line); (c) 
space-averaged squared vertical current (J|); (d) spectral indices a sma u (thin blue line) and 
enlarge (thick yellow line); (e) cancellation exponents Kj z (green thin line) and Kh c (thick red 
line). Vertical grey dot-dashed lines indicate the time of the M and X flares. The red dashed 
vertical line indicates the transition time t*. The gray area is the stationary interval used 
for the statistical study. 
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Fig. 6.— Magnetic power spectra estimated on February 12th at 22:24UT (circles) and on 
February 14th at 03:24UT (diamonds). Power law fits and the corresponding spectral indices 
are indicated. 
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H Jz , 2011 , Feb 14 , 03:24 UT 



Fig. 7.— The signed measure //j(Z) calculated for the vertical current density J z on February 
14th at 03:24UT. Color scale is arbitrary. 
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Fig. 8.— The scaling of the cancellation function x(l) calculated (a) for the vertical current 
density J z and (b) for the current hclicity h c , on February 12th at 22:24UT (circles) and 
on February 14th at 03:24UT (diamonds). Power-law fits are also indicated, along with the 
corresponding cancellation exponents. 
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t[h] 

Fig. 9.— Top panel: the time evolution of B rms in the time interval including three M and 
one X class flares. Vertical bars represent the instant of major flares eruptions. Time is given 
in hours, starting at the beginning of the whole observation. Bottom panel: time derivative 
dB rms /dt , highlighting the sharp changes occurred at the major flares. 
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Fig. 10.— Left panels: the correlation coefficients Pp{t) and /?s(r) at various time lags r, 
for the indicated paris of variables. The dot-dashed vertical line indicates the time of the 
maximum Spearman correlation. Right panels: the corresponding scatter plots at the time 
lag of the Spearman correlation peaks. 
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